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ABSTRACT 

The density profiles of dark halos are usually inferred from the rotation curves of disk galaxies based 
on the assumption that the gas is a good tracer of the gravitational potential of the galaxies. Some 
authors have suggested that magnetic pinching forces could alter significantly the rotation curves of 
spiral galaxies. In contrast to other studies which have concentrated in the vertical structure of the disk, 
here we focus on the problem of magnetic confinement in the radial direction to bound the magnetic 
effects on the H I rotation curves. It is shown that azimuthal magnetic fields can hardly speed up the H I 
disk of galaxies as a whole. In fact, based on virial constraints we show that the contribution of galactic 
magnetic fields to the rotation curves cannot be larger than ~ 10 km s _1 at the outermost point of Hi 
detection, if the galaxies did not contain dark matter at all, and up to 20 km s" 1 in the conventional 
dark halo scenario. The procedure to estimate the maximum effect of magnetic fields is general and 
applicable to any particular galaxy disk. The inclusion of the surface terms, namely the intergalactic 
(thermal, magnetic or ram) pressure, does not change our conclusions. Other problems related with the 
magnetic alternative to dark halos are highlighted. The relevance of magnetic fields in the cuspy problem 
of dark halos is also discussed. 

Subject headings: galaxies: halos — galaxies: kinematics and dynamics — galaxies: magnetic fields — 
galaxies: spiral 



1. INTRODUCTION 

The rotation curves of spiral galaxies, especially dwarf 
galaxies and low surface brightness galaxies, are used to 
derive the profiles of the dark halos which are very useful 
as tests for the standard cold dark matter scenario (e.g., 
de Blok & Bosma 2002). The usual assumption is that the 
dynamics of the neutral atomic gas is a good tracer of the 
radial gravitational force. It can be shown that although 
spiral waves and non-circular motions may contribute to 
produce substructure in the rotation curves (e.g., de Blok, 
Bosma & McGaugh 2003), they do not alter the rotation 
curves significantly for the galaxies selected in these in- 
vestigations. Galactic magnetic fields could also affect the 
gas dynamics in spiral galaxies (Piddington 1964; Peratt 
1986; Battaner et al. 1992; Sanchez-Salcedo 1997; Ben- 
jamin 2000; Beck 2002, 2003). 

The random small-scale component of the galactic mag- 
netic fields acts as a pressure, giving support to the disk 
and, therefore, leading to a rotation (slightly) slower than 
the gravitational circular speed. By contrast, some con- 
figurations of the large-scale magnetic field can give rise 
to a faster circular velocity. Nelson (1988) and Battaner 
et al. (1992) proposed the so-called "magnetic alternative" 
to dark matter in which the contribution of the magnetic 
pinch could be large enough to explain the observed ro- 
tation curves of spiral galaxies without the necessity of 
dark matter. This is an interesting possibility because 
some universal properties such as the disk-halo conspir- 
acy or the Hi-dark matter constant relation, first noticed 



by Bosma (1978), would have a natural explanation in this 
scenario (Sanchez-Salcedo 1996a), as well as the truncation 
of stellar disks (Battaner, Florido, Jimenez- Vicente 2002). 
However, after the publication of the paper by Battaner 
et al. (1992), a plethora of possible difficulties faced by 
the model was pointed out by different authors (Cudde- 
ford & Binney 1993; Jokipii & Levy 1993; Persic & Salucci 
1993; Vallee 1994; Katz 1994; Pfenniger, Combes, & Mar- 
tinet 1994; Sanchez-Salcedo 1996b, 1997). In an extensive 
review, Battaner & Florido (2000) argue that all the short- 
comings can be overcome and, therefore, magnetism can 
explain the way in which galaxies rotate, without the help 
of dark matter. 

A certain MHD configuration has the ability of produc- 
ing a faster rotation of the whole plasma in the system, 
i.e. magnetic confinement, when the net contribution of 
the magnetic fields to the virial theorem is negative (in the 
conventional notation). In an isolated system, it is well 
known that the magnetic field contribution to the virial 
theorem, applied to a large enough volume, is positive or 
zero, but never negative, reflecting the net expansive ten- 
dency of magnetic fields (e.g., Shafranov 1966). Jokipii 
& Levy (1993) used this argument to conclude that in the 
presence of magnetic fields, galaxies would need even more 
dark matter than in the unmagnetized case. One could 
argue that, since the virial theorem uses global (volume- 
integrated) variables, it is possible that magnetism has 
a "centrifugal" action in the vertical direction but cen- 
tripetal in the radial direction (radial confinement) and, 
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therefore, the expansive tendency of magnetic fields man- 
ifests only in the vertical direction. We advance that this 
claim is very wrong. An alternative way out is to assume 
that galaxies are not at isolation but embedded in the ubiq- 
uitous (magnetized) intcrgalactic medium which would be 
responsible for the required radial confinement. 

Motivated by these ideas and given the importance of a 
correct interpretation of the rotation curves for Cosmology, 
we will try to put bounds on the effects of magnetic fields 
on the H I rotation curves. However, we do not restrict our- 
selves to discuss the viability of explaining rotation curves 
without dark halos; the orbits of Galactic globular clusters, 
the satellite galaxies and possibly the dynamics of remote 
stars in the outer Galactic halo, which arc not affected 
by magnetic fields, strongly argue in favor of a dark halo 
around the Milky Way (e.g. Freeman 1997). The main 
issue is to see whether galactic magnetic fields can alter 
appreciably the H I rotation curves, especially in external 
galaxies for which the dark matter content is almost en- 
tirely derived from the rotation curve of the gas. Only if 
the answer is negative, we can be confident that the pro- 
files of dark halos, the dark matter content in galaxies as 
well as other interesting empirical correlations found be- 
tween, for instance, halo parameters and luminosity are 
reliable. 

Our investigation will be focus on the radial equilibrium 
configuration of magnetized rotating disks. For reasons 
that should become clear later, the analysis of the radial 
equation of motion puts stringent bounds on the maximum 
magnetic effects to the rotation curves; this being free of 
other complicating assumptions that appear in studies of 
the vertical equilibrium configuration such as the exter- 
nal pressure, the vertical dependence of the magnetic field 
within the height of the H I thin disk, or the flattening of 
the dark halo, which are all very uncertain. 

The paper is organized as follows. In §2 wc describe the 
ingredients of a simplified model in which magnetic fields 
could give rise to a faster rotation in the outer parts of H I 
disks. In §3 the maximum effect on the rotation curves 
is studied. Other difficulties inherent to the magnetic al- 
ternative are discussed in §4. The case of magnetic fields 
producing a slower rotation and their role as a possible 
remedy to the cuspy problem of dark halos is left to §5 
and the final remarks are given in §6. 

2. INGREDIENTS OF THE MAGNETIC MODEL AND BASIC 

THEORY 

In this Section, we describe some basic concepts useful 
for understanding the remainder of the paper. The issue 
of the problem of radial confinement will be addressed in 
§3- 

2.1. The One- dimensional Case 

In presenting the role of the magnetic tension in the 
H I rotation curves, the equations are simplified by assum- 
ing an equilibrium, steady-state, axisymmetric disk of a 
slightly ionized gas with a purely azimuthal large-scale 
magnetic field (e.g., Battaner et al. 1992; Sanchez- 
Salcedo 1997). The radial component of the equation of 
motion of the gas at z — in cylindrical coordinates reads 

V * {R) - R dR + -pdR + 8nRp dR ' (1) 



where p{R) is the gas density at z = 0, v^R) the cir- 
cular speed of the gas, &(R) the gravitational potential 
created by all masses in the system, and P(R) is the gas 
pressure consisting of the thermal, turbulent, convective, 
and cosmic ray components plus also the magnetic pres- 
sure arising from the random magnetic field component, 
but docs not include the magnetic pressure arising from 
the large-scale magnetic field. All these variables should 
be taken at the midplane of the disk. In equation (1) it was 
assumed 3> vr,v z , and thus we kept only the dominant 
terms associated with v^. In terms of the total pressure 
Ptot(R) = P(R) + B|(i?)/87r, it can be rewritten as: 
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The second term of the RHS is the radial pressure force 
that acts inwards or outwards depending on whether P tot 
increases outwards or inwards. The third term is the mag- 
netic tensile force which always imposes a net inward force. 

We should remind that ordered galactic-scale magnetic 
fields arc not expected to be a general feature in all disk 
galaxies, especially in dwarf galaxies. Although observa- 
tions of magnetic fields in these galaxies are scarce, they 
do not indicate magnetic fields with a global symmetry, 
but a much more chaotic galactic-scale field than in nor- 
mal spirals (e.g., Chyzy et al. 2000). Therefore, models 
based on a regular azimuthal magnetic field may not be 
an adequate representation for all of them. In fact, indirect 
observations suggest that regular magnetic fields may be 
ignorable, at least in some dwarf galaxies. The dominance 
of spherical holes in the Hi of the nearby dwarf galaxy 
IC 2574 has lead Walter & Brinks (1999) to conclude that 
large-scale magnetic fields are dynamically unimportant in 
the outer parts of this galaxy. 

As it is shown in Appendix A, P(R) is always a de- 
creasing function of R in approximately-isothermal, ex- 
ponential galactic disks, even if the disk were embedded 
in an ambient medium with high pressure. Therefore, 
the asymmetric drift always produces a force pointing 
outwards and hence it can never contribute to enhance 
the circular speed of the gas. Moreover, the asymmet- 
ric drift is almost negligible in rotating galactic disks ex- 
cept in very low-mass galaxies such as DDO 50 (Carignan, 
Beaulieu & Freeman 1990; de Blok & Bosma 2002; Bu- 
reau & Carignan 2002). Thus we will also ignore it, hence 
dPtot/dR w (l/8Tr)dB^/dR, and the force balance equa- 
tion is simplified to 



v 2 JR) = v 2 c (R) 



1 d{R?B\ 
8nRp dR 



(3) 



where v 2 = Rd^/dR is the gravitational circular velocity. 
Note that the magnetic pressure arising from the large- 
scale field is kept in our analysis because it cannot be ne- 
glected in general. The role of the magnetic pressure term 
associated with the random (small-scale) magnetic field 
component is discussed in detail in §5. It should be also 
clear that the vertical gradient of P is not negligible when 
considering the vertical configuration of the disk, since it is 
the main agent responsible for the vertical support against 
gravity. 
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From Eq. (3) we see that for a purely azimuthal mag- 
netic field decaying radially not faster than R~ x , the mag- 
netic tension force produces radial confinement. Such con- 
figuration is usually called a linear or Z-pinch. The re- 
sulting inward magnetic force, the pinch, accounts for the 
required centripetal force responsible for an anomalous 
faster rotation. Our aim is to constrain the maximum 
effect of the Z-pinch in galactic disks from very general 
grounds. It is worth mentioning that the Z-pinch configu- 
ration is not used in tokamaks as interchange and other in- 
stabilities quickly destroy the magnetic configuration (e.g., 
Roberts 1967). Restricted to purely azimuthal fields, the 
force-free magnetic configuration occurs for a field profile 
B(j, = Bq(Ro/R) 2 , where Bq is the magnetic strength at 
some radius i?o- For that profile, the outward magnetic 
pressure force is in balance with the inward tensile force 
at any radius. This profile is, however, rather unphysical 
because it is singular at R = 0. 

Observations of galactic magnetic fields suggest that the 
radial component of the field, Br, is comparable to B^ at 
least in the optical disk and hence neglecting Br may not 
be a good approximation. One could argue, however, that 
the magnetic pitch angle becomes very small in the out- 
ermost parts of galactic disks, where measurements of the 
field orientation are scarce. Since Br is not well-suited 
to produce radial confinement, the most optimistic and 
simplest situation, in order to maximize the magnetic ef- 
fects on the rotation curves, occurs when the field is taken 
as purely azimuthal. As we arc interested in putting up- 
per limits on the role of magnetic fields we will assume 
throughout this paper that Br = in the outer disk to 
give the model as much leeway as possible, as it is assumed 
in the magnetic alternative. 

An important contribution to P is the turbulent pres- 
sure. For a certain component of the interstellar medium, 
the averaged turbulent pressure is given by pv 2 /2, where Vt 
is the one-dimensional turbulent velocity dispersion. For 
the Hi components, Vt ~ 6-10 km s _1 (e.g., Boulares & 
Cox 1990) and is nearly constant with galactocentric ra- 
dius and also from galaxy to galaxy (Sellwood & Balbus 
1999, and references therein). The magnetic field is said 
to be in equipartition with the turbulent motions when 
the gas-to-magnetic field pressure ratio is 1, or equiva- 
lently, when the turbulent velocity dispersion coincides 
with the Alfven speed. In the literature, one often finds 
the statement that dynamo action in a thin disk saturates 
at equipartition field strengths (e.g., Ruzmaikin, Sokoloff 
& Shukurov 1988; Shu 1992; Poezd, Shukurov & Sokoloff 
1993; Beck et al. 1996; Kulsrud 1999). Therefore, it is 
interesting to estimate those values of v$ — v c for which 
the large-scale magnetic field would be in superequiparti- 
tion; this is done in the next. Denoting by va the Alfven 
velocity associated with the regular field, we can write 
v 2 - v 2 = {^R P )- 1 d{R 2 Bl)/dR = v\ + {R/8np)dBl/dR. 
In the likely case that the strength of the large-scale mag- 
netic field does not increase with galactocentric radius it 
implies v\ — v 2 < v\. Hence, for changes in the rota- 
tion velocity of v 2 — v 2 > v 2 , the regular magnetic field 
should be above equipartition with the turbulent motions. 
The above condition is equivalent to > \Jv\ + v 2 or, in 



terms of — v c : 

v^~v c > (v 2 + vl) 1 ' 2 - v c « v c (l + - v c = \ J. 

(4) 

For v c w 100 km s , variations of — v c > 0.5 km s 1 
already require magnetic fields over equipartition. Such a 
situation is not very appealing to most dynamo theorists. 
However, in order to present our results as generic as possi- 
ble and since there are some observations that could sug- 
gest magnetic energy densities above the turbulent ones 
in the disk of galaxies (e.g., Beck 2002, 2003) and in the 
Galactic center (e.g., Moss & Shukurov 2001), we will not 
employ equipartition when discussing the magnetic effects 
on the rotation curves. 

So far, the usual way to proceed was to integrate Eq. (1) 
out to a certain radius in order to derive the required radial 
profile for B${R) to explain the observed rotation curve vj,, 
given the volume density at z — and the gravitational 
circular velocity (including a contribution of the dark halo 
as in Sanchez-Salcedo (1997) or only the luminous com- 
ponents as in Battaner et al. (1992)). The question that 
arises is whether these local solutions can be matched with 
boundary conditions at infinity or not. This issue is very 
delicate and it will be proven very useful to put bounds on 
the influence of galactic magnetic fields in the gas dynam- 
ics at large scales. Before that, however, it is illustrative 
to describe how galaxies can reach a configuration with a 
faster rotation in order to have a complete picture of the 
underlying physics. This is briefly described in the next 
subsection. 

2.2. The Phase of Magnetic Contraction 

The assumption that the magnetic field is mainly az- 
imuthal immediately suggests that it commenced to grow 
after the galactic disk has formed. In the absence of 
magnetic tension, the initial galactic disk or protodisk 
should have a rotation speed given by the mass distri- 
bution w| = v 2 = Rd<&/dR. If the rotation curves are 
significantly affected by the magnetic field, one should ex- 
pect that, as the magnetic field is amplified no matter 
by what process, the disk initially in gravitational force 
balance, will experience a radial contraction caused by 
the growing magnetic tension, leading to a disk with a 
faster rotation due to angular momentum conservation. 
We will refer to this process as the magnetic contraction 
phase (MCP). If the mass loss by outflows is unimpor- 
tant, the MCP is the only way capable of spinning up the 
disk because the total angular momentum stored in the 
fields (oc J r x [E x B] dV) is negligible in both the initial 
protodisk and in the final configuration. Other processes, 
such as turbulent viscosity, are very inefficient to produce 
the desired effect. 

The conservation of angular momentum implies that, for 
a galaxy similar to the Milky Way with a dark halo and a 
rotation speed of 220 km s _1 at 30 kpc, an increase of 10 
km s -1 in the velocity speed requires a radial contraction 
of 1.3 kpc at that radius. The contraction must be huge 
in the magnetic alternative; see §4.1. 

3. GALACTIC DISKS WITH FASTER ROTATION. THE 
PROBLEM OF THE RADIAL CONFINEMENT 
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3.1. Isolated Galaxies 

In the previous section it was pointed out that there is 
a pinching effect only if B^ does not decrease faster than 
1/R. However, this decay is untenable at large R in real 
astrophysical systems as it will be shown in the next and, 
therefore, it is inevitable that beyond a certain radius, the 
magnetic tension produces a force pointing outwards (see 
Fig. 1 for a visualization). In fact, the slowest radial mag- 
netic decay occurs in the ideal case of having strict cylin- 
drical geometry with currents, J z , flowing from z = — oo 
to z — oo (otherwise, if currents are confined to a certain 
region, B must decrease with radius at least as fast as a 
dipolc field for large r: B < K/r 3 , where K is a constant). 
Under strict cylindrical geometry 1 , the condition of closed 

current loops is equivalent to imposing J Q fl J Z R' dR' — > 
at large R (hereafter, this will be our definition of isolated 
system in a cylindrical system). From Ampere's law, this 
requirement implies that B<p < C/R at large R, being C a 
constant. So there must exist a certain radius denoted by 
i?ko P at which the magnetic field starts to decay as 1/R or 
faster. R^ cp is therefore the radius at which the rotation 
speed of the gas coincides with the local Kcplerian velocity. 
If i?Hi is the radius of the last point of H I detection, we 
require i?kc P > Rm in order to have an anomalous faster 
rotation along the observed disk. 

The system will be in a configuration of equilibrium at 
R > i?kop only if the inward gravitational force is higher 
than the outward magnetic force 

d* 1 <K#Bj) 

dR 8-KRp dR w 
If this condition is not fulfilled the disk will expand ra- 
dially due to the magnetic forces. This condition can be 
written as 

1/2 



r>2 n2 

- n 4ccp- D 0,kcp 



8tt 



pRv 2 c dR 



-Rkc 



(6) 



where B^^ ep is the magnetic field strength at i?kop- Since 
a dependence B$ w C/R at large R is unphysical (see 
above), an upper limit to -S^.kcp can be established from 
Eq. (6): 



Bi 



,kcp 



< 



8tt 

^kep 



p(R, O)Rvl(R) dR, 



(7) 



nature of magnetic fields is expansive or null (but never 
confining) in both the radial and vertical directions (Ap- 
pendix B). More specifically, for toroidal magnetic fields 
and strict cylindrical geometry, the total rotational energy 
of an isolated system should be the same with or without 
magnetic fields, keeping the mass distribution fixed (so 
that <J>(i?); see Appendix B for further details). By impos- 
ing that constraint, a similar upper limit for i?ke P is found 
(Appendix C), except by a factor (X p b), between 1-2, that 
takes into account possible three-dimensional effects. 

It is convenient to consider first the extreme and hypo- 
thetical case that spiral galaxies do not host dark mat- 
ter halos. For a galaxy with a mass in gas and stars M, 
$(R) w -GM/R in the outer parts and 



f>2 



<Bi 



8ttA 



pB 



GM 

-n-kop 



p(R,0)dR. (8) 



The assumption that &(R) ~ —GM/R is satisfactory for 
many high surface density galaxies (e.g., van Albada et 
al. 1985) but it is not necessarily fulfilled in some low sur- 
face brightness galaxies (LSB) and dwarfs for which the 
gas surface density is not exponential and its contribution 
to the potential may be comparable or even higher than 
that from the stellar component at radii ~ Rm for cer- 
tain values of the adopted stellar mass-to-luminosity ratio. 
However, we are confident that our conclusions would not 
change for these galaxies even adopting the most unfavor- 
able assumptions. 

In the optimistic and rather unlikely situation that B^ 
were constant with R out to i?ke P , i.e., B$(R) = B cr at 
R < -Rkop, the characteristic variation — v% would be 
v \ ~ v c ~ B^KATip). By using Eq. (8), it is possible to 
estimate the maximum relative variation of the circular 
speed of the gas at a radius R < Rk cp : 



R 

Rkcr 



_Rg 

-Rkc 



cxp 



R — -Rkc 

R e 



(9) 

where we have adopted a radial exponential decay for 
p(R,0) with scale R K (e.g., Bland-Hawthorn, Freeman & 



Quinn 1997). Note that (vi 



;)/«c is not sensitive to 



Pext because even though the density at the midplane in- 
creases with P C xt, equation (9) does not depend on the 
absolute value of the density but only on R g . The main 
source of uncertainty stems from the value of i?kop, but the 
most optimistic situation occurs when i?ke P ~ Rm, so we 
will assume for the discussion that both radii are identical, 
~ 30 kpc for common spiral galaxies. In order to estimate 
v \ ~ v c from Eq. (9), we will assume that the mass of gas 
in the disk beyond Rm and within \z\ < h\, being h\ the 
characteristic scale height of the gas disk (see Appendix 
C for further details), is 0.5 percent the total mass of the 
disk. This seems to be the case of dwarf galaxies (van 
den Bosch et al. 2001), but it is probably a too generous 
value for normal spiral galaxies. The latter assumption 
corresponds to a radial scale length for the gas, R g , of 
~ 5.5 kpc. Putting these values into Eq. (9), we obtain 
( v l ~ v c)/ v c = 0.36A pB at R = 30 kpc. If we take for 

1 In such a situation, the magnetic field is always purely azimuthal, consistent with the assumptions given in §2 and, hence, the essential physics 
of the mechanism giving rise to greater azimuthal velocities of the H I disk, is entirely contained within the one-dimensional radial models. 



where Rb^ is the radius where B^ w 0. In practice we will 
take Rb^ — *• oo to estimate the maximum contribution. 
The maximum value for -B^kep corresponds to a situation 
in which there is balance at R > Rk cp between gravita- 
tional and magnetic forces and thus the material outside 
the i?kop circle has no rotation at all. 

The upper limit derived above has an easy interpreta- 
tion. In principle, one could achieve very large differences 
Vc/, — v c if one would make the density arbitrarily small than 
the Alfven speed would be arbitrarily large. However, ac- 
cording to Eq. (5), this is not possible globally because the 
magnetic field has to be held in by the radial weight of the 
plasma in the radial gravitational field since the net effect 
of magnetic fields is never confining. 

It is customary and natural to discuss this kind of con- 
straints using the virial theorem which says that the net 
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der S of height L, can be written as: 



2T+2 / r<i\ •_'/'< I • IV • .V-l, , 



2 

(10) 

where T is the kinetic energy of the gas, W its gravi- 
tational energy in the gravitational field created by all 
masses in the system, P the gas pressure as defined in 
§2.1, Ps the gas pressure at the boundary surface, r is 
a coordinate vector, I = d 2 I/dt 2 with I the trace of the 
moment of inertia tensor and Ai C yi is given by 



M 



cyl 



reference values v c = 100 km s _1 at 30 kpc, the maxi- 
mum change in the rotation velocity, Av = — v c , ranges 
between 15 to 30 km s _1 for \ p b = 1 to 2, respectively. 

This value for Av should be regarded as a generous es- 
timate for two reasons. First, the azimuthal field is not 
expected to be approximately constant while the volume 
density of gas at z — varies by a factor of <~ 50-100 
(Diplas & Savage 1991). Secondly the adopted equilib- 
rium configuration, namely v^(R) = at any R > i?ke P , 
cannot be reached in practice, starting from a rotating 
disk in equilibrium with the gravitational force, due to the 
conservation of angular momentum. In fact, the magnetic 
field begins to grow in the disk as discussed in §2.2, creat- 
ing a magnetic force directed outward beyond i?kep- This 
region will undergo a continuous outward drift, while the 
inner disk (R < Rk ep ) must be subjected to the required 
contraction phase, generating a gap in the surface density 
of the disk. In order to quantify this effect, we assume for 
simplicity that the outer disk (i? > i?kcp) suffers an ap- 
proximately self-similar expansion with factor ip > 1, i.e., 
p(R) — (po/'0 2 ) exp(— R/ipR g ) at R > V'-Rkcp, in the final 
equilibrium configuration. From the conservation of angu- 
lar momentum, it can be shown that the right-hand side of 
Eq. (9) is now a factor ip 2 /(ip — 1) smaller (Appendix C). 
Hence, the optimized situation occurs for ip — 2 and then 
( v l - v c)/ v c < O.lAps at R = 30 kpc. This implies that 
Avtj, < 10 km s _1 at the periphery of Hi disks in normal 
spiral galaxies under the assumption that they contained 
only stars and the observed gas. Therefore, it is not pos- 
sible to banish dark halos by magnetic fields because one 
would require much larger values of Av^. 

For disk galaxies with a dark halo, as postulated in the 
conventional scenario, it is possible to follow the same pro- 
cedure to estimate the maximum effect of magnetic fields 
when v c =cte, instead of v 2 = GM/R (see Appendix C). 
Similar arguments suggest that the expected local, and 
probably transient, variations in the rotation speed larger 
than <~ 20 km s _1 are very unlikely even in models includ- 
ing dark halos (see also §2.2 and Sanchez-Salcedo 1997). 

In the next subsections, it will be shown that these up- 
per limits for Av^ also apply even when galaxies are not 
isolated systems but are embedded in an ambient medium, 
thereby making these estimates very robust. 

3.2. Non-Isolated Galaxies 

It could be argued that the model described by Eq. (1) 
is too simplistic because other terms not considered in that 
equation, as a radial flow or other components of the mag- 
netic field, might be important, regarding the radial con- 
finement, beyond -Rhi 2 - In fact, unlike the radial compo- 
nent of the field, the z-component may contribute to the 
radial confinement of the disk and will be included in our 
analysis. In addition, the hydrodynamical pressure by a 
radial inflow due to intergalactic accretion of matter may 
also play a role. Since systems embedded in an ambient 
medium usually present sharp boundaries, it is convenient 
to use an integral form of the equations as that given by 
the virial theorem. 

The scalar virial theorem for fields with Br = and 
strict cylindrical symmetry (i.e., d/dz = d/d<j) = 0) ap- 
plied to a finite cylindrical volume V , bounded by a cylin- 

2 Within the Rhi circle, observations suggest that other terms are unimportant for those galaxies under consideration. 
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(e.g., Fiege & Pudritz 2000 and Appendix B). Through- 
out this paper the subscript S refers to quantities valuated 
at the radial surface of our cylindrical volume with radius 
i?s, and the subscript 'cyl' stands for the cylindrical case. 
This form for the virial theorem is valid for any arbitrary 
helical field {Br = 0) provided cylindrical geometry; no 
other additional assumption is required. The surface in- 
tegral terms including the z = ±L/2 surfaces are taken 
into account in Eq. (11). The integration of the magnetic 
terms can be found in Appendix B. For any further details 
in the derivation, we refer the reader to the paper by Fiege 
& Pudritz (2000). 

We can identify two new vias of radial confinement: 
by ram pressure [last term of LHS of Eq. (10)] or by 
an external magnetic field. The ambient pressure (Ps), 
on the other hand, only produces vertical confinement as 
discussed in Appendix A. If these new terms are able to 
produce radial confinement then the constraint given by 
Eq. (7) could be somehow relaxed. The possible role of 
these terms will be discussed in two separate sections. 

3.2.1. Ram Pressure 

So far, we have neglected the terms involving vr, be- 
cause H I observations show that vr is very small compared 
to v,), for the galaxies under consideration. Beyond the 
edge of the observed Hi disk, the value of vr is uncertain. 
We consider now the surface integral of the momentum 
stress in Eq. (10), which may have a negative sign for the 
case of strong mass accretion into the disk. In fact, it can 
be shown that even though I is expected to be negative, 
the combination — 1/2 — J s p (v ■ r)(v ■ dS) could be still 
negative, showing the confining effect of the ram pressure. 
(For example, in the one-dimensional case of a supersonic 
flow that slows down to rest by crashing perpendicularly 
to a rigid plane wall and forms a strong adiabatic shock, 
-J/2 -J sP (vr)(v dS) = -(2/3) J s p(v r)(v ■ dS).) 

Since spherical accretion cannot confine a disk, the in- 
flow should be extremely collimated parallel to the galac- 
tic plane, i.e. vr ^> v z , for which there is no evidence. 
Moreover, this collimation is very unlikely to occur if the 
disk forms from a sphere of gas. Gas cools and parcels of 
gas acquire circular orbits, conserving approximately an- 
gular momentum, when they collapse to settle into a disk, 
leading to vr — > (see Rognvaldsson 1999 for results based 
on numerical simulations). Furthermore, in the hypothet- 
ical case that such inflow were possible by some unknown 
mechanism, the accretion speed should be so high that an 
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accretion shock would be unavoidable, transforming ram 
pressure into thermal pressure and probably generating a 
hydraulic jump in the edge of the disk and other observable 
features such as Ha emission. The hydraulic jump will di- 
vert the flow around the disk decreasing dramatically its 
confinement ability. 

3.2.2. Confinement by the Intergalactic Magnetic Field 

A coherent external (intergalactic) magnetic field could 
in principle contribute to the confinement under certain 
circunstances. It is important to distinguish here between 
the intracluster and the wider intergalactic field. As we 
are dealing with galaxies in the field, we are interested in 
the amplitude of the latter one. At present, setting obser- 
vational constraints to the strength and topology of the 
intergalactic magnetic fields is a challeging task (see Ryu, 
Kang & Biermann 1988; Kronberg ct al. 2001; and Widrow 
2002 for a review) and, therefore, they are very uncertain. 
Coherent magnetic fields on scales much greater than the 
typical sizes of galaxies could be as large as <~ 0.1 fiG in 
the optimistic situation that extra amplification through 
batteries, dynamos or other routes of magnetization were 
operative. 

We assume the ideal situation that initially the pregalac- 
tic magnetic field is uniform and fills the protogalaxy mak- 
ing a finite angle a with the rotation axis of the galaxy. 
As the disk is forming by accretion of gas, the magnetic 
field lines are advected with the flow. It is well-known in 
studies of the collapse of interstellar clouds that, under 
flux-freezing conditions, this compression produces that 
the total magnetic contribution to the virial theorem is 
positive (see, e.g., Chapt. 24 of Shu 1992 and Chapt. 11 
of Mestel 1999) and, therefore, the confinement by the ex- 
ternal field is impeded. In order to have confinement, the 
disk must be able to get rid of the magnetic flux. 

In the strict cylindrical symmetry, the external magnetic 
field must be vertical and then Eq. (11) with = 
tells us that magnetic confinement is only possible if the 
volume-averaged vertical component of the magnetic field 
within V is strictly smaller than the magnetic field at the 
boundary. Since vr is directed inward and Ovr/OR < 
in the MCP, \B Z \ increases at every point of the disk due 
to the inward advection (see also Reyes-Ruiz & Stepin- 
ski 1996; Moss & Shukurov 2001); under those conditions 
confinement is not viable. This problem could be allevi- 
ated if magnetic field lines were expelled from the disk by 
ambipolar and turbulent diffusion. However, a thorough 
numerical study by Howard & Kulsrud (1997) shows that 
\B Z \ is also increased by ambipolar shearing terms to a 
value considerably larger than its initial value. In the lack 
of a mechanism able to reduce \B Z \ within disks, pregalac- 
tic magnetic fields do not help to confine galactic disks. 
This example clearly illustrates why self-confinement is so 
rare in nature. 

Models based on the confinement by anchored fields to 
the galactic nucleus will present a similar problem, as \B Z \ 
is expected to decline outwards. 

4. THE MAGNETIC ALTERNATIVE 

Battaner & Florido (1995, 2000) and Battaner, Florido 
& Jimenez- Vicente (2002) pointed out that: (1) the mag- 
netic alternative to dark halos does not suffer from the 



problem of excessive flaring of the H I disk and (2) the ex- 
istence of truncated stellar disks is explained satisfactorily 
in the magnetic alternative. Putting aside the problem 
of radial confinement discussed in the previous sections, 
it is stimulating to see whether alternative scenarios to 
conventional ideas have the potential of explaining other 
phenomena in a natural way. If this occurs, alternative 
models could attrack much more attention. For instance, 
some empirical correlations of galactic warps have led to 
Castro-Rodriguez et al. (2002) to suggest tentatively that 
outer rotation curves could not be caused by the pres- 
ence of a massive halo. In the next sections we will show 
that the flaring problem remains unsolved and that the 
existence of truncated stellar disks cannot be considered a 
prediction of the magnetic alternative. 

4.1. The MCP and Implications for the Truncation of 
Stellar Disks 

If galaxies contained only stars and the observed gas, the 
initial protodisk of mass Mi should follow approximately a 
Keplerian rotation, — GMJR, at large R, and evolves 
to the present-day exponential disk of mass Mf and radial 
scale length i?d, due to magnetic contraction. The peak 
velocity of the final disk is v max = 0.62y/ GM { /R d (e.g., 
Binney & Tremaine 1987). For simplicity, we assume that 
beyond 2.27 i?d the rotation curve of a present galaxy is flat 
with circular velocity v max - During the MCP, a ring of gas 
located initially at a radius i?; must squeeze up to a final 
radius Rf, such as v max Rf — ^/GMLJli from angular mo- 
mentum conservation. Therefore, the initial radius is given 
by the simple expression R, = 0.38(M f /M i )(.Rf/- R d)- R f- 
For a normal spiral galaxy like the Milky Way, with a 
nearly flat rotation curve up to the last measured point, 
typically at <~ 30 kpc, and with i?d ~ 4 kpc, the MCP 
should be able to squeeze the disk at least from <~ 85 kpc 
to ~ 30 kpc in radius, provided the mass loss in outflows 
is negligible, corresponding to a typical radial velocity of 3 
km s _1 for one Hubble time. According to §3, this contrac- 
tion is an unfeasible process due to the global expansive 
nature of magnetic fields. 

The level of contraction may be significantly reduced if 
a substantial fraction of the mass is lost in winds following 
early epoch of star bursts, i.e. Mf < M\. Recently, Zhao 
(2002) has estimated that less than 4% of the total mass 
of the Galaxy has been lost during its lifetime. Even as- 
suming an amount of mass loss 3 times larger than this 
value, the magnitude of the contraction of the disk is still 
huge, from 72 kpc to 30 kpc, again for i?d = 4 kpc. 

The MCP represents a serious threat to the suggestion 
by Battaner, Florido, Jimenez- Vicente (2002) that, since 
stars do not feel the galactic magnetic field, the escape of 
stars formed beyond a certain galactocentric radius could 
explain the observed truncation of stellar disks. The cal- 
culations by Battaner et al. (2002) are rooted in a scenario 
that does not consider the MCP, with the gas disks being 
in anomalous rotation from the very beginning. If one 
tries to include the MCP in the model, the stellar dis- 
tribution function depends on different assumptions that 
concern the details of the MCP, such as the radial velocity 
of the contraction, or the relationship between star for- 
mation rate, surface density and magnetic field strength. 
Consequently, the escape radius depends on time. If the 
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disk was contracted from ~ 80 to 30 kpc in radius, it could 
be even more natural to construct a model in which the 
stellar surface density decays with R slowly enough, with 
no stellar truncation radius. 

Even assuming that stars in the outer disk were not 
formed until the gas disk stopped its contraction, which 
is in contrast with observations that show old stars in the 
outer disk (e.g., Ferguson & Johnson 2001), the planar ve- 
locity dispersions of the stars should increase dramatically 
with R and this is not observed (e.g., Pfenniger et al. 1994; 
Vallee 1994). 

4.2. The Flaring Problem 

In the past years, much of the challenge involved in mak- 
ing the magnetic alternative viable concerned the exces- 
sive flaring of the disk, which is caused by the vertical 
magnetic pressure (Cuddeford & Binney 1993). In a two- 
dimensional analysis, Battaner & Florido (1995) predict a 
scale height of 3-3.5 kpc at a radius of 15 kpc and of 7 kpc 
at 30 kpc for the H I disk of M31. These values are still too 
large by a factor of 3-5 if comparing with the observed H I 
scale height observed in the Galaxy by Diplas & Savage 
(1991) and Burton (1992). 

This discrepancy in the scale height is the least of our 
worries. An inspection of the radial profile of the magnetic 
strength (fig. 3 of Battaner & Florido 1995) reveals that 
there is something deeply wrong in their calculations. As 
shown in §2.1, the magnetic field cannot decrease faster 
than 1/R in order to have a disk with a faster rotation. 
Therefore, a magnetic field of 5 /iG at 15 kpc implies a 
field strength larger than 2.5 /iG at 30 kpc. The magnetic 
field estimated by Battaner & Florido (1995, 2000) decays 
faster than 1/R. 

If the dark matter contribution to the rotation curve 
is null, the required magnetic fields in the outer parts 
of typical spiral galaxies are so intense (~ 5-10 /iG) that 
the problem of excessive flaring is unresolved. The flaring 
problem is exacerbated when the magnetic pressure of the 
small-scale magnetic field is included because it also con- 
tributes to the vertical expansion of the disk. This pressure 
is observed to be comparable to the magnetic pressure of 
the regular field at least in the optical disk of the Galaxy. 
On the other hand, the excessive flaring is also aggravated 
if the magnetic pitch angle is not zero. We conclude that 
the flaring problem by itself is sufficient to reject the mag- 
netic alternative. 

5. GALACTIC DISKS WITH SLOWER ROTATION: THE 
CUSPY HALO PROBLEM 

In this Section, we will discuss the role of magnetic fields 
as a remedy for the cuspy problem of dark halos. It is 
interesting to see whether the magnetic pressure of the 
random component, as well as certain configurations of 
a global-scale magnetic field (not necessarily azimuthal) 
could help to explain the discrepancy between the profiles 
of dark halos obtained in cold dark matter (CDM) simula- 
tions, which becomes very cuspy towards the center (e.g., 
Navarro, Frenk & White 1996, 1997), and observations of 
the halos of dwarf galaxies and LSB, which are better fit- 
ted with a core-dominated halo (see de Blok, Bosma & 
McGaugh 2003, and references therein). One possibility is 
that galaxies are indeed embedded in cuspy halos follow- 



ing the profile found by Navarro, Frenk and White (1996; 
hereafter NFW) but the combination of asymmetric drifts 
and magnetic effects wipe out completely the inner cusp, 
leading to rotation curves that mimic a halo with a core. 

The magnetic pressure of the random component gives 
some support to the disk, reducing the rotation speed. 
Since it is expected to be in cquipartition with the gas 
turbulent pressure, corrections due to the magnetic pres- 
sure should be of the order of the asymmetric drift, ~ 2-3 
km s _1 (e.g., de Blok & Bosma 2002), which are smaller 
or comparable to observational uncertainties, typically 4-6 
km s _1 . This contribution, however, could be larger in a 
certain radial portion of the disk if the magnetic pressure 
decays faster than vt ■ Writting the Alfven speed associated 
to the random field as v\ = g(R)v%, with the turbulent ve- 
locity dispersion v t constant with R, and dg/dR < 0, the 
rotation curve changes according to the relation: 

9 9 /„n v ? f dhip d\nq\ 

*=* + * R H(idi + jdi)- (12) 

Following de Blok & Bosma (2002) , we define R ln as the ra- 
dius of the innermost sampled point of the rotation curve, 
and assume that at R > Ri n , g{R) = {R m /R) n for sim- 
plicity. Doing so, we are imposing that at i?; n , «a ~ v t . 
The second term of the RHS of Eq. (12) is of the order 
of the asymmetric drift (comparable to observational un- 
certainties), whereas the last term is — nv\ R n \/2R n . As 
the observed rotation curve differs from the NFW pro- 
file along a significant range of galactocentric radius, say 
at R 3> R[ n i=a 500 pc, then < n < 1, implying that 
nv%R?j2R n < vl/2 ~ 50 km 2 s~ 2 (for v t = 10 km s" 1 ). 
Consequently, this effect would be only significant very 
close to the galactic center where v c < 10 km s _1 , a re- 
gion with large observational uncertainties. Summing up, 
the contribution from the total pressure to the circular ve- 
locity is less than 5-6 km s _1 . This result agrees with 
the fact that CO, Ha and Hi rotation curves reveal the 
same kinematics for the central regions of some galaxies 
(e.g., Bolatto et al. 2002), which suggests that magnetic 
fields do not contribute significantly to the support of the 
gas disks. In particular, the CO rotation curves are ex- 
pected to trace the motion of the highest density regions, 
including molecular clouds, which are thought to be hardly 
affected by the galactic magnetic fields. 

In conclusion, magnetic effects are not sufficient to solve 
the cuspy halo problem satisfactorily. 

6. FINAL REMARKS 

It has been recognized for decades that the galactic mag- 
netic field provides significant vertical support to the in- 
terstellar medium. The role of the magnetic tension in the 
vertical hydrostatic configuration was considered by Cox 
(1988) and Boulares & Cox (1990), whereas the role of 
the tension in the plane of the disk was studied by Nel- 
son (1988), Battaner et al. (1992) and Benjamin (2000) 
[see also Piddington 1964]. The idea that large-scale mag- 
netic fields can explain the rotation curves of spirals with- 
out invoking dark matter has been revived by Battaner & 
Florido (2000) and Battaner et al. (2002). 

All the previous work aimed at giving bounds to mag- 
netic effects on gas dynamics are based on the flaring of the 
disk by the vertical magnetic gradients. Here we propose 
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that the radial analysis provides stringent constraints, be- 
ing free of other complicated assumptions such as the un- 
known external pressure. We have shown that magnetic 
fields have a net expansive (or null, but never confining) ef- 
fect not only in the vertical direction but also in the radial 
direction. Without dark matter halos, the contribution of 
the magnetic field to the rotation curves at the periph- 
ery of spiral galaxies is not expected to be larger than 10 
km s -1 , which is insufficient to banish dark halos by mag- 
netic fields. This result is valid even though an external 
pressure, intergalactic magnetic fields or ram pressure are 
taken into account. In addition, the huge contraction of 
the disk required in the magnetic alternative would spoil a 
simple magnetic explanation for the observed truncation of 
stellar disks. Moreover, the problem of excessive flaring in 
the outer disk by the magnetic pressure is unresolved and 
may be aggravated if the random component of the mag- 
netic field is included, or if the field is not purely toroidal. 
By contrast, there is observational evidence that the dwarf, 
dark-matter dominated galaxy IC 2574 does not possess a 
dynamically-important galactic-scale magnetic field. Our 
first conclusion is that the magnetic alternative is not vi- 
able. 



Even though there is no doubt about the existence of 
dark matter halos around galaxies, it is desirable to con- 
strain the magnetic effects on the H I rotation curves for a 
correct interpretation of the observations. It is suggested 
that magnetic fields can contribute to produce wiggles in 
the Hi rotation curves, which could be likely associated 
with non-axisymmetric waves in the disk, as may be the 
case of NGC 1560. However, the magnetic tension is un- 
likely to produce substructure in the rotation curve of mas- 
sive spiral galaxies above 20 km s _1 . In Sanchez-Salccdo 
(1997), the detailed shape of the rotation curve of this 
dwarf galaxy was fitted reasonably well combining a mag- 
netic field of strength ~ 1 /j,G plus an isothermal dark 
halo. Finally, we have considered magnetic fields as a pos- 
sible solution to the cuspy problem of dark halos. We 
conclude that they are insufficient to explain the discrep- 
ancy between the observations and the predictions of CDM 
simulations. 

We thank J. Canto and A. Hidalgo for discussions, 
and the anonymous referees and the Editor for critical 
comments. F.J.S.S. acknowledges financial support from 
CONACYT project 2002-C40366. 



APPENDIX 

A. WHY AN EXTERNAL PRESSURE IS UNIMPORTANT TO STUDY THE RADIAL CONFIGURATION 

Consider a three-dimensional axisymmetric gas disk with surface density £(P) supported mainly by rotation and 
embedded in an ambient medium of pressure P ex t- The transition between the gas disk and the external medium occurs 
at a vertical height zs(R)- From the vertical hydrostatic equilibrium, the radial gradient of the total (turbulent plus 
magnetic) gas pressure Ptot at the midplane of the disk is given by 



dPtot(fl,0) 
dR 



d_ 

dR 



-Pext + 



L 



z s (R) 



<9$ , 

Pir dz 

oz 



d_ 

dR 



z s (R) 



<9$ , 

P d; dZ 



(Al) 



where $(R,z) is the gravitational potential. The pressure exerted by the ambient medium squeezes the disk along z, 
increasing P tot (P, 0) but barely dP to t{R, 0)/dR, and thereby making the radial confinement extremely difficult. This can 
be visualized in the following manner. Ignoring for a moment magnetic fields, an external pressure will produce some 
radial confinement if it is possible to reach a configuration with dP/dR\ > just by increasing P ext : 

' rzs(R) q$ 
Jo P ^ dZ 



_d_ 

dR 



> 0. 



(A2) 



For a gas density decaying along z with a scale height h, for instance p oc exp(— z 2 /2h 2 ) (for \z\ < zg), Eq. (A2) is 
satisfied if the scale height increases with galactocentric radius at least as fast as h ~ ~ irE -1 , where I is an exponent 
that varies from 2-3 for a gravitational potential with a flat and a Keplerian circular speed, respectively. However, this 
requirement cannot be fulfilled in approximately isothermal galactic disks with the surface density decreasing outwards, 
because the maximum variation of h occurs in the limit of large external pressure and goes as h ~ R l . 

There are other physical difficulties in the boundaries of the disks if they were confined by P cx t- From dimensional 
grounds, it is expected that the gas pressure P(R, 0), comprising mainly thermal and turbulent contributions as defined in 
§2.1, should be greater than the external pressure, as a consequence of turbulent motions and mixing by Rayleigh- Taylor 
type instabilities in the disk. The total pressure at z = is 

Ptot {R, 0) = P(R, 0) + ^ > P oxt + Q • (A3) 

In the periphery of the disk, pressure confinement requires P to t(P, 0) w P ex t, thus Eq. (A3) implies P|(P, 0)/87r <C P ext < 
P(R, 0), i.e. well below equipartition and, therefore, they are not sufficient to affect appreciably the rotation curves, 
according to our discussion in §2.1. 



B. DERIVATION OF THE MAGNETIC TERMS IN THE VIRIAL THEOREM 

Here we discuss the magnetic terms that appears in the virial theorem. The same derivation can be found in Fiege & 
Pudritz (2000). Consider a distribution of a electrically conducting plasma with velocity v in a gravitational field $(P) 
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created by all masses in the system. The scalar virial theorem applied to a finite volume V, bounded by a surface S, can 
be written as: 



,d 2 ' 

o (v ■ r)(v ■ dS) = —- 

is 

where T is the kinetic energy of the plasma within V, W the gravitational term given by W = — J v pr ■ V$ dV and 



2T + U + M cyl +W- J p(vr)(vdS) = ^, (Bl) 



IV i 

U = 3 [ PdV- [ Pr-dS, (B2) 
Jv Js 

M = lj v BidV + ±-j s (r.B)(B. dS) - JL ^ B 2 r ■ .5, (B3) 

/ pr 2 dV 7 (B4) 



I : 



'V 

where r is a coordinate vector (e.g., Shu 1992). 

In the case of having cylindrical symmetry (i.e., d/dz = d/dcj) — 0) and for fields with no radial component, Br = 0, a 
more simple expression for M cy i, which denotes M under cylindrical symmetry, can be derived when the virial theorem 
is applied to a finite cylindrical volume; this is done in the following. Since Br — 0, the second term in the RHS of 
Eq. (B3) contributes only at z = ±L/2 surfaces, where L is the vertical length of the volume V. The component of the 
field perpendicular to these surfaces is B z , so that 

- J^r • B)(B .dS) = - I 2irRB 2 z dR = - | B 2 dV. (B5) 

The last term of Eq. (B3) comprises an integral along the vertical surface of the cylinder I\ and a contribution of the cap 
surfaces, I2, given by 

1 



and 



h = -^(Bl,s + Bl s )V, (B6) 

T f Rs 1 f 

I 2 = / 2ixRB 2 dR = / B 2 dV. (B7) 

8?r J a 8?r J v 



Combining the three contributions the final expression for M C yi reads as 

^cyi = ^(J v B l ^ (Bi s + B 2 tS ) Vj . (B8) 

Similar manipulations allow us to rewrite U as U = 2 J v PdV — 2PsV. 

For a cold galactic disk, the thermal energy is negligible compared to the kinetic energy. If, in addition, there is no 
radial accretion flows and the disk can be assumed to be in a steady state, the virial theorem is reduced to 

2T + M cy i + W = 0. (B9) 

Assuming that galaxies are isolated systems, B z ,sV and B^^V go to zero at large R implying that always .Mcyi > 0, 
reflecting the impossibility of magnetic confinement in the radial direction. In fact, for a given distribution of mass, W is 
fixed and if M cy i is positive, the total kinetic energy should be smaller than in the unmagnetized case in order to satisfy 
the virial constraint. This does not exclude the possibility that some parts of the disk could present a faster rotation 
(larger kinetic energy) as long as rotation is slower than that given by gravity in other parts of the disk. It might be that 
we are observing just the part of faster rotation and the material with a slower rotation lies outside the last point of H 1 
detection, if it is fully ionized by the metagalactic UV background (see Fig. 1 for a visualization). For the case of purely 
azimuthal fields M cy i = 0, so 2T + W = 0. Therefore, given a distribution of mass, all the configurations in equilibrium 
must have the same T regardless of magnetic fields. Note that M cy i — does not imply necessarily that the magnetic 
configuration is force-free. In the particular case of a force-free magnetic configuration, it can be shown that Ai cy i — at 
any arbitrary volume and hence magnetic fields are neither expansive nor confining. 



C. CONSTRAINING THE AZIMUTHAL MAGNETIC FIELD STRENGTH IN ISOLATED GALAXIES 

Consider an isolated disk galaxy comprising of gas and stars with a given total surface density E(i?). The gas is moving 
in pure circular orbits (i.e. v z = vr = 0) and threaded by a purely azimuthal magnetic field. The distribution of mass 
(stars+gas+dark matter) will create a gravitational potential at z = 0, <I>o(-R)- The gravitational circular velocity (not 
the observed rotation velocity) is defined as v 2 = Rd&o/dR. We will assume for simplicity that along the thickness of 
the Hi disk the rotation curve is constant, i.e. dv<f ) /dz\n = 0. The disk is supposed to be rotating with an anomalous 
rotation speed v^R) > v c (R) out to the radius i?kc P - 

We split up the rotational energy, T, in two parts T = Tq + AT, where To is the kinetic energy as if the disk, as it 
is, were rotating in force balance with the gravitational potential <&o without magnetic fields (keeping all the rest of the 
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variables of the disk S(i?) and so on unchanged). More specifically, T between two arbitrary radii, say R\ and i?2, and 
within \z\ < hi, where hi is the characteristic scale height of the disk, is given by: 



Ri 




Ri i-R 2 r 




= T 


+ 27T / / 


Ri 




Ri JRi Jo 



pR(vl-v 2 c )dz dR, 



(CI) 



where p{R,z). According to our discussion in Appendix B, M. cy \ = for pure azimuthal fields, and thus the total 
rotational energy of the disk must be the same with or without magnetic fields. Therefore, AT — between R = and 
oo, implying 

rRkc 
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Since we have assumed that v^R) > v c (R) at R < i?kc P , this equation implies that v^R) < v c (R) outside the i?kc P circle. 
Using Eq. (1), the LHS of the equation above reads 
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where -B^kep is the field at i?kcp an d A^ is a dimensionless coefficient less than 1 defined by the relation J 1 B 2 (Ryr Cp , z) dz = 
XshiB^R^cpjO). Analogously, the RHS of Eq. (C2) can be written as: 

poo phi poo 

/ / P R(v 2 -v 2 )dzdR = X p h 1 p(R,0)R(v 2 -v 2 )dR, (C4) 

•> Rkcp JO J Rkcp 

where A p is a dimensionless parameter less than 1. Thus, the value of £?0,kcp turns out to be 
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The factor \ p s = A p /As takes into account possible three-dimensional effects. For a disk with a constant scale height 
with radius, one would expect X p /\b ~ 1, whereas for a disk that flares, it could be a little bit larger, depending on the 
value adopted for hi. As hi — > 0, X p /Xb goes to 1, whereas for a value of order of the scale height, ~ 400 pc, X p /Xb ~ 2 
could be considered as an optimistically high value. 

The maximum magnetic effect on the rotation velocity within the i?kep circle occurs when -B^.kcp takes its maximum 
value; this happens when v<f,(R) = for R > R^ cp 
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B cr is an upper limit to the magnetic field strength at i?ko P - 

In the standard dark matter scenario and for galaxies with a flat gravitational circular speed, v c (R) =constant, the 
maximum relative contribution of the magnetic field, assuming generously that B^{R) = B cr at R < i?kep (see §3.1), is: 



vl - v 2 



2A 



pB 



& 

-Rkcp 



1 



Rkcp 



exp 



R — R 



kcp 



(C7) 



where an exponential law with scalelength R g was adopted for p{R, 0) beyond i?kc P - 

In the other extreme case that no contribution to the gravitational potential comes from a dark component, v 2 (R) ~ 
GMj R at large R, where M is the mass of the luminous matter in the disk. The corresponding maximum relative variation 
is given by 
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If the disk is initially, in the lack of magnetic fields, rotating with velocity v c , the growing magnetic tension during 
the MCP will produce a contraction of the inner disk and a radial drift of the outer disk. Let us assume that the disk 
beyond i?kcp expands in a self-similar fashion and achieves a final density profile pf(R,0) = (po/ip 2 ) exp(— R/ipR g ), with 
tp > 1 the scale factor. Tthe conservation of angular momentum implies that v c / > (R) = ip~ 1 v c (R/ip), where v c (R/tp) is the 
circular velocity at R/ip. In the conventional dark halo scenario, v c ^constant so that v<p(R) = v c /ip. Substituting these 
relations into Eq. (C5), we get 

_ 8^ 2 - 1)A^ C 2 



B 



,kcp 



iP 2 R 2 



kcp 



p { (R,0)RdR. 



'A -Rko; 



Therefore, the relative effect on the rotation curves starting with a disk in gravitational force balance is 



2A p s 



V - 1 \ f R, 

~ 02 



1 



_Rg 

^?k< 



exp 



R — R 



kcp 



Rcr 



(C9) 



(C10) 



\ Rkcp / \ ± ^kep 

Similar manipulations allow us to calculate the same relation but assuming that galaxies do not contain massive dark 
halos. The corresponding calculation will give 



ip-1 



R 



Rcr 



R-Rke 
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Fig. CI. — Schematic diagram showing the effect of the magnetic forces in a vertical cross-section of the gas disk in the magnetic alternative 
(panel a). R s i zc is the radius of the gas disk. Beyond RhIj the outermost radius with Hi detection, the gas could be fully ionized. The big 
arrows indicate the net effect of the magnetic fields onto the gas disk. Beyond Rkepi the magnetic forces point outwards in all directions. 
The expected rotation velocity versus radius is shown in panel b (solid line), together with the gravitational circular rotation speed given by 
baryonic matter (dashed line). In both panels R has the same scale. 



